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Abstract 

We calculate the decay rate ratios for OZI allowed decays of T(5S') to two B mesons 
by using the decay amplitudes which incorporate the wave function of the T(5S') 
state. We obtain the results that the branching ratio of the T{5S) decay to B*Bg 
is much larger than the branching ratio to BsB^ or BgB*, in good agreement with 
recent experimental results of CLEO and BELLE. This agreement with the exper- 
imental results is made possible since the nodes of the T{5S) radial wave function 
induce the nodes of the decay amplitude. We find that the results for the T{5S) 

decays to B^B'^' or B'J'B^ ' pairs are dependent on the parameter values used for 
the potential between heavy quarks. 



PACS codes: 12.39.-x, 12.39.Pn, 13.25.-k, 13.25.Gv 
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1 Introduction 



The CLEO Collaboration recently observed the Bg mesons in the e'^e annihilation 
at the T{5S) resonance [1]. They established that Bg meson production proceeds 
dominantly through the creation of B*Bg pairs. They found a{e~^e~ — B*Bg) = 
[0.1llo:o3(stat) ± 0.02(syst)] nb and set the following limits [1]: 



a(e'^e B^B^ 



a(e+e- ^ BtB,] 



< 0.16 



[a(e+e- ^ bX) + ^(e+e" ^ B^ ^ ^ _ 



a(e+e- ^ B*B^) 



The BELLE Collaboration collected data at the T{5S) resonance and obtained the 
following ratio at the T{5S) energy [2]: 



aje+e- ^ B^X) 



93+^ ± 1 % . (3) 



The above results of CLEO and BELLE show that the branching ratio of the T(5S') 

^ 

decay to B*Bg is much larger than others ones to Br'B\ ' pairs, and in particular 
that the ratio of \a{e^e~ bX^ + o-(e+e" BlBs)\ to a-(e+e" B^l) is much 
smaller than the value which is expected from the spin countings. These results imply 
that there exists a delicate mechanism in the decays of T{5S) — > BB. 

The unitarized quark model [3] could predict similar results for the above ratios. 
These ratios were also calculated in Ref. [4] by using the decay amplitudes which 
take into account the bound state effect through the wave function of the T{5S) 
state. Even though these theoretical calculations showed the nature of the above 
experimental results, they were performed before the masses of the B mesons, espe- 
cially Bg and B*, are measured. Therefore, it is desirable to perform a theoretical 
calculation by using the experimentally established mass values of all B mesons in 
order to obtain the results which are more reliable quantitatively. It is the purpose 
of this paper. We use the same method for the calculation as that used in Ref. [4] 
with the experimentally measured values of B meson masses. We perform the numer- 
ical calculations for two sets of the parameter values for the potential between heavy 
quarks, in order to see the dependence of the results on the parameter values. 

The CLEO Collaboration also measured the cross sections of the T(5S') decays to 
B*B*, BB*, and BB, and found that the decay to B*B* is dominant with (74± 15)% 
of the total B rate [5]. Here, B = Bu or Bd, and BB signifies both BB and B*B. 
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They summarized the measured branching ratios as follows [1, 5, 6]: 

B{B*B*) = 43.6 ± 8.3 ± 7.2 % , 
B{BB* + B*B) = 14.3 ± 5.3 ± 2.7 % , 

B{BB) < 13.8 % , (4) 

B{Bi*^B''*^) = 16.8 ±2.6+11% (using yields) , 

B{Bi*^B^*'^) = 24.6 ± 2.9j:^^3° % (using yields) . 

The above branching ratios show that B{B*B*) is about two or three times of 

B{Bi*^B^*^). In this paper we also compare the results of our calculations with the 
above experimental results. 

This paper is organized as follows. In section 2 we briefly summarize the mecha- 
nism of light quark pair creation and the formula of the decay amplitude. In section 
3 we explain the procedures and results of our calculations. The last section is con- 
clusion. 



2 Formula for Decay Amplitude of T{bS) BB 

The Cornell group studied the effect of OZI allowed decay channels [7]. They proposed 
that the following interaction hamiltonian is responsible for the decay as well as the 
binding of quark-antiquark bound states [7]. 

k=i ^ 

where 

V{r) = -- + -,. (6) 

In (5) Pfc(x) = ^/'^(x)|Afc^/'(x) are the color densities of quark fields. This model 
corresponds to the vector coupling since (5) is the leading term of the vector coupling 
hamiltonian in the nonrelativistic expansion. 

The decay rate of T(5S') with spin s (which is 1) and mass M to 5 mesons with 
spin si and S2 is given by the formula [4] 

T = ^^E,E,C{ss,S2)\A,s{P)\\ (7) 

where P is the magnitude of the center of mass momentum of the mesons and Ei 
and E2 are their energies. The spin counting factor C{ssiS2) are for s = 1 given be 
1/3, 4/3, and 7/3, for the cases Si = S2 = 0, Si = S2 = 1 + Si = 1 S2 = 0, and 
■Si = S2 = 1, respectively [7]. 
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Fi gure 3: I^q(^P) for the potential (A). The sohd hne is for the Bs meson pairs 
with TJig = 0.55 GeV, and the dotted hne is for the Bu^ or B'^^ meson pairs with 
= = 0.33 GeV. 

The decay amphtude A^siP) calculated from the diagrams in Fig. 1 is given by 
[7, 8, 9, 10] 

A,s{P) = 4vr/,/iJP) , /, = {-f (m.^^^^y ; q = u,d,s, (8) 

V vr / V rrib + niq/ 

where for the vector coupling of (5) the momentum dependent function lli^{P) is 
given by [7, 8, 9, 10] 

iLiP) = dtmRnd^m^) (9) 

with 

m = [te-^\it^-l)e-^'/'^erfi^)]+W4-te-'\e-''/^ (10) 

where Rnhif^) is the radial wave function of the heavy quark-antiquark bound state, 
is the spherical Bessel function, and pb = nib/ i'm-b + mq). We note that lliiP) 

incorporates the radial wave function i?„L(r) of the T(5S') state. The first and second 
terms in (10) come from the linear and Coulombic parts of (6), respectively. Following 
[7], we use the radial wave functions of the final cu and cd states with L = which are 
given by the ground state harmonic-oscillator wave functions (47r)2 (^)4exp(— /3r^). 
For the value of f3 in the calculations of (9) and (10), we use /3 = (2^)(|^)^ [7], 
where /i is the reduced mass rriqpb of the cu, cd or cs system. 
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Figure 5: The position of the final momentum P of the B meson for each mode 
in the T(5S') rest frame, (a) is for the potential (A) and (b) is for the potential 
(B). 
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T(5S) 


Bu 


Bd 


b: 


B*, 


Bs 


b: 


10865 


5279.15 


5279.53 


5325.1 


5325.1 


5366.3 


5412.8 



Table 1: Meson Masses (MeV) used in the calculation [12]. 



BuBu BuB^ 


BIB., 


BdBd BdB^ 


BdBd 


BsBs BgBg 


b:b. 


1.282 1.183 


1.075 


1.280 1.182 


1.075 


0.846 0.682 


0.462 



Table 2: Decay momenta in GeV for T(55') calculated with masses given in Table 1. 

3 Results of Calculation 

We calculate the decay rates for OZI allowed decays of T(55') to two B mesons by 
using the formula (7) with the decay amplitudes given in (8) and (9). For V{r) in (6) 
we use the following two parameterizations in our calculation: 

(A) Eichten et al. [7]: k=0.517, a=2.12 GeV"\ 771^=5.17 GeV. 

(B) Hagiwara et al. [11]: k=0.47, a = 1/^019 GeV-\ mb=4.75 GeV. 

For the masses of u, d and s quarks, we use m„ ^=0.33 GeV and ms=0.55 GeV. 

The radial wave functions of T(5S') obtained by solving the Schrdinger equation 
with the parameter sets (A) and (B) are presented in Fig. 2. The function IIq{P) in 
(9) obtained with the above (A) and (B) sets of the parameter values are presented 
in Figs. 3 and 4. 

In order to obtain the decay amplitude for each final state of B meson pair from 
(8) and (9), we need to know the final momentum P of the B meson in the T(5S') 
rest frame. We use the meson masses presented in Table 1, which give the value P for 
each mode as written in Table 2. In Fig. 5 we marked the position of the value P for 
each mode, from which we can get the magnitude of /|o(-P) each mode of B meson 
pair. Then, we can calculate the decay rate F from the formula (7) for each mode of 
the B meson pair. Using the values of F obtained in this way, we calculate the ratios 
of the decay rates. We present our results of this calculation in Table 3. When we 
compare our results presented in Table 3 with the experimental results in (4), for the 

Bd'^^d 

pairs we should use two times the values of Table 3 since in (4) B can be 

(*) — (*) 

either neutral or charged [5, 6]. In our calculation the results for the Bu B^ pairs 

are the same as those for the S^*^!?^*"* pairs since we use the same value 777.^^^=0.33 
GeV for the u and d quark masses. 

We find in Table 3 that the results for the Bi*"^ B^*^ pairs are in good agreement 
with the experimental results in both cases of (A) and (B) sets of the parameter 
values. Especially, we note that the result for the ratio [cr(e+e~ — > -Bs-B*) + cr(e+e" — > 
B*Bs)]/(T{e^e~ — B*Bg) is consistent with the experimentally obtained very small 
value. The reason why we could obtain such a consistent result is that we incorporate 
the radial wave function of the T(5S') state which has nodes. The nodes of the radial 
wave function induce the nodes of the decay amplitude as we can see in Figs. 3 and 
4. 

We find that the results for the T(55') decays to Bu^ B^*^ or B'^^ B^^^ pairs are 
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Potential 




BsBg + BgBl 


BgBs 


B^B^a 


BdB^ + BdB*^ 


BdBd 


(A) 


L + C 


1 


0.007 


0.065 


3.159 


0.298 


0.036 


(A) 


L only 


1 


0.010 


0.059 


3.142 


0.312 


0.030 


(B) 


L + C 


1 


0.034 


0.141 


0.923 


0.294 


0.646 


(B) 


L only 


1 


0.026 


0.130 


0.953 


0.242 


0.579 



Table 3: Decay rate ratios of T(5S') BB modes. Potentials (A) and (B) are given 
in the text. "L + C" and "L only" mean that we used "both linear and Coulomb 
parts" and "only linear part" respectively for the calculation of (10). We always used 
both linear and Coulomb parts for the calculation of the wave function of T(55'). 



sensitive to the parameter values used for the potential between heavy quarks. For 
example, the result with the set (A) for the ratio [B{BlE*^) + B{B*B*a)]/B{B*Bl) 
is rather larger than the experimental result and the result with the set (B) for this 
ratio is within the error bar of the experimental result. Such a situation happens 
since the node positions of the decay amplitude are sensitive to the parameter values 
of the potential in the range of large final B meson momentum. 

4 Conclusion 

The CLEO and BELLE Collaborations recently performed experiments at the T(5S') 
resonance. Their results for the ratio [a{e~^e~ —>■ B^B^) + a{e^e~ B*Bs)]/a{e^e~ — *• 
B*Bl) at the T(5^) energy is much smaller than the value expected from the spin 
countings. This delicate phenomenon was predicted by the unitarized quark model 
[3] and by the decay amplitudes which incorporate the wave function of the T(55') 
state [4, 7]. However, these theoretical studies were performed before the masses of 
the B mesons, especially B^ and B*, are known. We use the same method for the 
calculation as that in Ref. [4] with the experimentally measured values of B meson 
masses. We perform the numerical calculations for two sets of the parameter values 
for the potential between heavy quarks, in order to see the dependence of the results 
on the parameter values used. 

Our results for the ratios for the Bi*^B^*^ pairs are in good agreement in both 
cases of (A) and (B) sets of the parameter values of the potential between heavy 
quarks. It is especially important that the result of our calculation agrees well with 
the delicate experimental result that the branching ratio to BgB^ or B^B* pair is 
much smaller than that to B*Bg pair. Obtaining such results by our calculation was 
possible since we take into consideration the radial wave function of the T(55') state 
which has nodes, and the nodes of the radial wave function induce the nodes of the 
decay amplitude of the T(5S') BB reaction. However, our results for the T(55') 

decays to Bi*^B^*^ or ^J^^J^ pairs are sensitive to the parameter values used for the 
potential between heavy quarks since the node positions of the decay amplitude are 
sensitive to the parameter values in the range of large final B meson momentum. 
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